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Abstract 

The Yang-Mills functional integral is studied in an axial variant of 't Hooft's 
maximal Abelian gauge. In this gauge GauB' law can be completely resolved 



^ ■ resulting in a description in terms of unconstrained variables. Compared to 

previous work along this line starting with work of Goldstone and Jackiw 
one ends up here with half as many integration variables, besides a field 
living in the Cartan subgroup of the gauge group and in D-l dimension. 
The latter is of particular relevance for the infrared behaviour of the theory. 
Keeping only this variable we calculate the Wilson loop and find an area 
law. 
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I. INTRODUCTION 



It is a common belief, that the fundamental interactions are described by gauge theories. 
This is, in particular, true for strong interactions, which are assumed to be described 
by QCD. This theory has been tested in the high energy regime, where perturbation 
theory is applicable due to asymptotic freedom. On the other hand, low energy hadron 
physics requires a non-perturbative treatment of QCD. This regime, which is alternatively 
related to the confinement problem, is much less understood. Perturbative calculations 
indicate, that the confinement phenomenon is due to the non-abelian nature of Yang- 
Mills theory. Furthermore, one-loop calculations show, that the perturbative Yang-Mills 
vacuum is unstable Q and various models of the Yang-Mills vacuum have been designed 
as, for example, the Copenhagen vacuum the instanton liquid model , which extends 
the instanton gas picture the dual superconductor || or the stochastic vacuum ||. 
The various models aim at different aspects of strong interaction, e.g. the instanton 
models seem to be suited to explain spontaneous breaking of chiral symmetry 0, while 
the Copenhagen vacuum, stochastic vacuum or the dual superconductor focus on the color 
confinement. 

A rigorous approach to strong coupling Yang-Mills theory is provided by lattice Monte- 
Carlo calculations ||, which have been developed to a high level of sophistication. This 
approach has given much insight into the nature of the Yang-Mills vacuum. The great 
successes of lattice calculations are in low-energy hadron physics (where confinement is 
perhaps not of much relevance) [[|. 

There has been much progress in controlling finite size effects. This concerns both the 
fermion description on lattice as well as improved lattice actions for gauge fields |10|| . 
However, a complete understanding of the Yang-Mills theory will probably not be provided 
by the lattice simulations alone but requires also analytic tools. For some applications 
of lattice QCD a separation of scales is required and input from perturbation theory is 
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needed (see e.g. [fLlp . Also the interpretation of the lattice results sometimes requires or 
at least is facilitated by modelling properties of the Yang-Mills-vacuum like correlation 



functions, which in turn are fed by the lattice calculations (see e.g. [O])- 



Several analytic approaches have been proposed to explore the non-perturbative features 
of strong coupling Yang-Mills theory. In particular, there is the strong coupling expansion 
itself for lattice QCD ||. In the strong coupling limit one finds Wilson's area law (or a 
linear rising quark potential), which is considered as a signal of confinement. However, 
the strong coupling limit yields the area law for almost any gauge group and in any 
dimension, but the finite convergence radius of the strong coupling expansion forbids the 
strong coupling result to be applied to the continuum limit. In fact the finite convergence 
radius is taken as an explanation to argue away confinement in cases where the property 
of confinement is undesirable, as e.g. in QED. 

Luscher at al. [|13|] have proposed a small volume expansion exploiting the fact, that 



perturbation theory becomes eventually applicable as the volume decreases [14 



Most analytic approaches to Yang-Mills theory are based on the Weyl gauge A = 0, 
where the GauB' law has to be inforced as a constraint to guarantee local gauge invariance 



17J . Violation of GauB' law generates color charges during the time evolution and, by 
this leaking of color, confinement is lost. This fact has recently been emphasized in 
reference [[L8[], where explicit projection on gauge invariant states has been performed in 
the construction of the path integral representation of the Yang- Mills transition amplitude. 
Not surprisingly projection onto gauge invariant states requires a compact integration 
measure (over the Haar measure of the gauge group), reminiscent to the lattice approach. 



In fact, the approach of ref. [18| can be obtained from the lattice formulation by taking 



the continuum limit in the spatial directions only. 

Several approaches have been advocated, which explicitly resolved the GauB' law con- 
straint by changing variables resulting in a description in terms of a reduced number of 
unconstraint variables. These approaches are based on the Schrodinger functional for- 



mulation of Yang-Mills theory ||17|| . Refs. | l9fl , f20|j use variants of unitary gauge, while 
in ref. J2J]] the Coulomb gauge VA = was used. Refs. |0|, |20| basically end up in a 
description in terms of gauge invariant variables (Further approaches along these lines are 



proposed in [22], [23 



Recently an alternative descriptions of Yang-Mills theory in terms of gauge invariant vari- 
ables constructed either from the magnetic |PJ] or electric J25] fields have been proposed. 

1, has 



In ref. [2^] the long wave length (strong coupling) limit of the formulation of ref. [| 
been studied, exploiting methods from the description of collective excitations of atomic 
nuclei. 

A crucial point in all analytical approaches to Yang-Mills theory is the gauge fixing, which 
in non-abelian gauge theories cannot be performed in a unique way due to the existence of 



Gribov copies [27]. 't Hooft has proposed a so-called maximum abelian gauge ||15|| , which 
partly eliminates the non-abelian components of the gauge field (abelian projection). In 
this gauge monopoles occur and it is believed, that their condensation yields confinement 
via the dual Meissner effect, see also |23. This approach has been further elaborated 
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16fl , where the first lattice calculations along this line have been performed. In fact, 



recent lattice calculations performed in the maximum abelian gauge show evidence for 
monopole condensation [^ . 



Recently, QCD on a spatial torus has been considered in the Weyl gauge A = in the 
canonical quantization approach [30]. Using a (partial) axial gauge, which, in fact, is a 
variant of a maximum abelian gauge, the resolution of GauB' law has been achieved by 
applying unitary gauge transformations, which rely on the quantum field operators. This 
results in a Schrodinger description in terms of unconstraint variables, where the resulting 
Hamiltonian is non-local. It is fair to say, that at the moment we have little experience 
with solving functional Schrodinger equations in quantum field theories. Furthermore, 
even in the functional (operator) approach matrix elements are given by D-l-dimensional 
functional integrals. Therefore it might be more convenient to use the (D-dimensional) 



functional integral representation from the very beginning. In the present paper I perform 
a resolution of Gauss' law analogous to ref. but in the functional integral represen- 



tation of Yang-Mills theory. It is the hope, that the functional integral formulation will 
facilitate in finding appropriate approximation schemes. Furthermore, the functional inte- 
gral approach provides more direct excess to the topological properties of the Yang-Mills 
vacuum and to numerical simultations, exploiting Monte Carlo techniques. 

The balance of the paper is as follows: In order to set notation and conventions, in 
section 2 the functional integral description of Yang-Mills theory is briefly reviewed and 
some relevant features are discussed. In sect. 3 and 4 we fix the gauge and resolve the 
GauB' law constraint. The Faddeev- Popov determinant is evaluated in sect. 5. In sect. 6 
the Wilson loop is evaluated thereby including only the dominant infrared unconstrained 
degrees of freedom. In sect. 7 the electric field variables are integrated out resulting in 
a theory in unconstrained degrees of freedom of the gauge potential. A short summary 
and some concluding remarks are given in sect. 8. Some calculations are relegated to 
Appendices. 

II. HAMILTON! AN FORMULATION OF GAUGE THEORIES 

We consider the gauge group G = SU (N) with anti-hermitian generators T a satisfying 
the commutation relation 

rpa rpb jabcrpc ^2 ^) 

where f abc are the structure constants. We choose the standard normalization 

tr (r a T b ) = --^5 ab . (2.2) 
Later on we will also make use of the generators in the adjoint representation defined by 

fa\ = _jabc / 23 X 
/ be 



which satisfy the same commutation relation ( |2.1| ). Throughout the paper we shall indi- 
cate the adjoint representation by the hat ,, A ". For a quantity \ a living in the gauge 
group we define the fundamental and adjoint representation, respectively, by 

X = X a T a , X = X a f a . (2.4) 

We also use the Cartan decomposition of the gauge group 

G = H®G/H, (2.5) 

where H = [/(l)^ -1 denotes the Cartan subgroup of the gauge group and G/H is the 
corresponding coset space. Furthermore, we use the indices with index zero, a , bo, . . . 
to denote a generator of the Cartan subgroup T a ° G H , j^o^fco = -while the 
indices a, b, . . . are reserved for generators of the coset space, G/H. Accordingly the gauge 
potential Ai(x) is decomposed as 

Ai = A™ + Af , (2.6) 

where Af = A^°T a ° is the gauge potential of the Cartan subgroup H and Af 1 lives in 
the coset space G/H. With respect to (color) charges of the Cartan subgroup H, Af is 
neutral, while Af 1 is charged. 

We also introduce the covariant derivative by 

Dp = d ll + A ll (2.7) 

and the field strength tensor 

F, v = [D^D u } = F« v T a 

F; v = d,A a u - d v A% + f ahc A\Al . (2.8) 
Under a gauge transformation Q(x) G SU(N) the gauge potential transforms as 

A^ — ► A% = tt (p^) = QA^ + n (fyfit) . (2.9) 
In the Hamilton formulation of Yang-Mills theory, which is based on the Weyl gauge 



A (x) = , 



;2.io) 



the dynamical variables are the spatial components of the gauge potential A a { . We shall 
use spatially periodic boundary conditions for the field variables 



A1 (x + L 



fce 



Ai{x) 



(2.11) 



where e k denotes 3-dimensional (spatial) unit vector, so that we consider Yang-Mills 
theory on a 3-dimensional torus. We have not yet specified the boundary condition in 
time direction. 

Let \C) denote an eigenstate of Ai(x), i.e. Ai(x)\C) = Ci(x)\C) where Cj(x) is a classical 
field function. The gauge invariant transition amplitude between static initial and final 
field configurations A4 (xq — 0, x) — C[ (x) and A4 (xq — T, x) — C" (x) is defined by |3l| , 



Z[C",C) = (C"\e- HT P\C) 



where 



H = Jd 3 x (^E«(x)E?(x) + ^ 2 Bt{x)B«{x)} 



(2.12) 



(2.13) 



is the Yang-Mills Hamiltonian with bare coupling constant g, electric field E%(x) = 
S/iSAl(x) and magnetic field B^(x) = ^e^jF^Ax). Furthermore P is the projector onto 
gauge invariant states 



G 



(2.14) 



Here is the vacuum angle |17j and the functional integration with respect to the Haar 
measure /i(O) of the gauge group extends over all time-independent gauge transformation 
Q n (x) with winding number n. For a gauge transformation Q (x) the winding number is 
defined by 

1 



24vr 2 



d xe ijk tr (RiRjR k 



R k = Qd k Q j . 



;2.i5) 



As usual we assume here that the gauge function fl (x) approaches a unique value floo for 
| of | — > oo, so that R 3 can be compactified to S 3 and n[Q] is a topological invariant. 

For many purposes it is sufficient to consider the partition function 

Z = J VCi(C\e- HT P\C) , (2.16) 

which can be easily reduced to the standard form 

Z = J2e~ EkT (2.17) 

k 

with Ek being the energy eigenvalues. Using the completeness of the eigenstates \k) of 
H {H\k) = Ek\k)) and P 2 = P it can be rewritten as 

Z = [ VC t £ ^ k (C)e- E ^* k (C) , (2.18) 

J k 

where 

* fc (C0 = (C\P\k) (2.19) 

are the gauge "invariant" energy eigenfunctionals, which under a gauge transformation 
Q n with winding number n transform as 

^ fc = e- ine ^ k {C) , (2.20) 



as is easily inferred from the explicit form of the projector ( 2.14 ). Assuming proper 
normalization of the energy eigenfunctionals ^(C), i.e. 

J Vd%{C)%{C) = S kl (2.21) 

eq. ( |2.18|) reduces to the standard form ( |2.17| ). 



In ref. it was explicitly shown that the gauge invariant partition function Q2.16j ) is 
given by the standard functional integral representation 



Z = J VA^gfe-^^^ ^ , (2.22) 



where 



is the usual Yang-Mills action and 



s ym [a\ ='j <^xf; v {x)f; v {x) (2.23) 



M\ = 3^ / d A xF; v Ff v (2.24) 

is the Pontryagin index with F* = \e^ UK \F K x being the dual field strength. The functional 
integration runs over all temporally periodic gauge field configurations (x = T) = 
■A-ii (xo = 0) and it is understood that the gauge fixing is included by the Faddeev-Popov 
method as indicated in fl2.22|) by S g f. 



At first sight one may wonder that eq. ( |2.22| ) reproduces the gauge invariant partition 



function ( 2.16| ) inspite of the missing Haar measure?]. However, as explicitly shown in 



the Haar measure arises from the Faddeev-Popov determinant. Similar investigations 
have been previously performed in ref. pTf . 



The equivalence proof between eqs. Q2.16|) and ( 2.22 ) |33] relies only on the gauge invari 



ance of the Hamiltonian and holds therefore also true when fermions are included. In this 
case the partition function is given by 

Z = J VqVq J T)A^^- SYM[A]+iev[A] , (2.25) 

where the fermion fields satisfy antiperiodic boundary conditions q (xq = T) = 
—q(xo = 0). For later convenience we rewrite the partition function as 

Z = J VqVqeI^Z YM [J] , (2.26) 

where 



1 In ref. IT^] it was claimed that the conventional functional integral representation ( |2.22| ) falls 
short of guaranteeing gauge invariance in the non-perturbative regime. 
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Z YM [J] = J VA^e- SYM[A]+ S J ^ +iBu[A] (2.27) 

is formally the partition function of a gauge field coupled to an external color current 
= qYJfiQ- Eq. (|2.27|) defines the Lagrange representation, which is fully covariant. 
For subsequent considerations it is more convenient to use the Hamilton formulation which 
arises from fl2.27| ) by linearizing the (F oi ) 2 term by means of an integration over the electric 
field variable Ef{x) which in view of eq. (|2.11|) has to satisfy spatially periodic boundary 
condition Ei(x + Lei) = E^x). Then the A field can be integrated out yielding the 
Gaufi' law constraint 

J VA expi J d A xAl(x)Y a (x) = 5 (T a ) , (2.28) 

where 

T{A, E) = diEi + [Ai, Ei) + J = [A, A] + Jo (2.29) 
is the generator of infinitesimal gauge transformations. Eq. (|2.27| ) then becomes the 



Hamilton functional integral representation of the partition function of Yang-Mills theory 
in the presence of an external source J M = (J°, J 1 ), which after continuing to Minkowski 
space and assuming = readsQ 



Z[C", C, J] = fv{A l1 E i )\{5U a {A,E))5{Y a {A 1 E))X{DetM a \x,) 

x xo 



ex p 1 3 / d * x 



E^A, - i (E?E? + BfB?) - AJ, 



(2.30) 



2 In fact in the derivation of the path integral representation ( 2.27 ) the Hamilton formulation 
(|2.30| ) arises in an intermediate step of the calculation. In the present case the Hamilton and 
Lagrange form are obviously completely equivalent. But in more general cases (e.g. in theo- 
ries with momentum dependent masses) the Hamilton form is obviously the more fundamental 
representation and the Lagrange form may even not exist. Furthermore the path integral deriva- 
tion shows, while the integral over the gauge field configuration has to be taken with periodic 
boundary conditions A4 (x° = T,x) = Ai (x° = 0,x), the integration over the electric fields is 
not constrained by any temporal boundary condition. 
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Here, V (Aj, Ei) denotes the (flat) functional integral measure over the gauge potential 
A® and the electric field Ef. Furthermore, f a (A, E) — is the gauge fixing constraint 
and DetAi ab (x,y), where xq = yo, is the Faddeev-Popov determinant. 

In the following two sections we will explicitly resolve the GauB law constraint S (T a ) 
and the gauge constraint 5 (f a ) in fl2.3U| ) leaving a functional integral over unconstrained, 
gauge fixed variables. 



III. GAUGE FIXING AND PARTIAL RESOLUTION OF GAUC LAW 



GauB' law ( |2.29| ) T a = has the generic form 

VE a = p a , p- = -[A h E % } a -J a ^ (3.1) 
where p a is the total color charge density. Applying GauB' integration theorem it follows 

dtE a = Q a , Q a = f d 3 xp a . (3.2) 



8M 

For periodic electric fields the electric flux through the surface of the box, § dY>E, vanishes. 
Consequently periodic boundary conditions to Ef(x) can only tolerate a vanishing total 
charge 

Q a = . (3.3) 

For the resolution of GauB' law a proper choice of gauge fixing is crucial. In the past a 
complete resolution of GauB' law has been achieved in the gauge e aik E® = for SU(2) 



in ref. [|l9|and an extension to SU(3) was considered in ref. [20]. There have been also 



attempts of a complete resolution of GauB' law in the Coulomb gauge [gl|] . The Coulomb 
gauge, which is singled out in QED by the absence of radiation of static charges has 
proved, however, to be inconvenient in non-Abelian gauge theories, in particular for an 
explicit resolution of GauB' law. In this respect axial type of gauges are much more 
convenient as was already realized in refs. ||34]| , [RH] and recently discussed in detail in 
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ref. ||30|| , where an explicit resolution of GauB' law in the canonical quantization approach 
has been performed. Below we will perform an analogous resolution of GauB' law in the 
functional integral approach. For this purpose it is convenient to choose the 3-axis as 
the preferred direction of the axial gauge and divide the GauB' law generator into parts 
parallel and perpenticular to the 3-axis 

T(x) = D 3 E 3 + T ± , r_L = D ± E ± + J . (3.4) 

If D 3 were regular the GauB' law T = could be easily resolved leading to an elimination 
of E 3 . Unfortunately, as we will explicitly see below, on the torus D 3 has always zero 
modes, independently of the used gauge. In fact, since D 3 transforms gauge covariantly its 
eigenvalues are independent of the gauge. Nevertheless we can exploit the gauge freedom 
to cast D 3 in as simple a form as possible. From this point of view the axial gauge A 3 = 
would be preferable. However, this gauge condition conflicts with the periodic boundary 
condition. This can be easily seen by considering the Polyakov line operator 

V 3 {x) = Pexp (<j> dx' 3 A 3 (x, x' 3 )^j , (3.5) 

where P denotes path ordering and the integration runs from a point[] x = (x, x 3 ) along 
the 3-axis to the point x = (x, x 3 + L). Due to the periodic boundary condition on A 3 the 
integration in ( |3.5| ) runs over a closed loop but nevertheless due to the path-ordering V 3 [x] 
depends on the starting point x. Under gauge transformation this quantity transforms as 

V 3 [x) — ► V^[x) = Q(x)V 3 {x)Q\x) (3.6) 

and one can obviously choose a gauge in which V 3 (x) is diagonal 

ff(x) = e as(x)L , a 3 {x) = a c 3 °T c ° . (3.7) 

However, it is impossible to gauge transform V 3 (x) to V 3 (x) = 10. 



3 Here and in the following we use the convention (x) = (x°,x) = (x,x 3 ) . 

This can be also easily seen in the lattice formulation. Starting at x 3 = one can bring 
the links U 3 {x) = exp (— aA 3 {x)) to the gauge U 3 {x) = 1 except for the last link terminating at 
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For the resolution of GauB' law it is convenient to follow ref. pO] and use the gauge 



A?(x) = 



or Al(x) = 



(3.8) 



This condition, of course, does not fix the gauge completely but allows still for arbitrary 
abelian gauge transformations u(x) G H . We will later make use of this freedom. Let us 
also mention, that the gauge transformation necessary to bring a given gauge field Ai(x) 



into the form (3.8) requires in general also topologically non-trivial gauge transformations. 



In the gauge ( |3.8| ) the operator D 3 is block diagonal with respect to the color neutral and 
charged components 



D 



ab 



Q §aobo o 



V 



Df 



(3.9) 



since f a ° b ° c = 0. Hence in this gauge the neutral part of the GauB' law generator simplifies 
to 



T a °(x) = d 3 E%° +T a f 



(3.10) 



On the space of periodic functions £ n (^) = e !IJ " X3 , u n = -j 1 the operator d 3 has a zero 
eigenvalue (n = 0) corresponding to a ^-independent eigenfunction. For simplicity of 
notation we have set here L = L 3 . The corresponding projection of E£ onto this zero 
mode 

L 



e 3 (x 



L 



dx 3 E^ (x, x 3 ) 



(3.11) 



does not enter T n (x) (|3.4j) and is hence not restricted by GauB law. 



Since the eigenfunctions of d 3 belonging to zero and non-zero eigenvalues are orthogonal 
in the Hilbert space of periodic functions, the neutral part of the GauB' law constraint 



£3 = L, which cannot be gauged away due to the periodic boundary condition. 
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separates in the two independent constraints corresponding to the subspaces of the zero 
and non-zero eigenvalues. Defining 

L 

-y± = jJd X3 ri , r7 = ri- 7 ±, (3.12) 


we have 

5(T n ) = 5 (d 3 E'« + I*) 5 ( 7± ) = 5 (d 3 E' 3 n + I*) 5 ( 7± ) , 

where 

E' 3 = E 3 -e 3 (3.14) 

lives entirely in the subspace of eigenfunctions with non-zero eigenvalues of d 3 . The 
constraint of the first 5-function can be easily resolved. Defining by d 3 the operator 
resulting from d 3 when the zero eigenvalue is removed, we obtain (with d 3 E' 3 n = d' 3 E' 3 ) 

i ^ + r" 1 ) = ^(£r + ir» 1 ). (3.15) 



(3.13) 



Hence the neutral part of GauB law eliminates the variable E' 3 . In addition we now 
exploit the residual invariance under Abelian gauge transformations to remove also the 
corresponding conjugate field variable 

A' 3 = A 3 -a 3 , a 3 = jJdx 3 A n 3 (3.16) 

by imposing the gauge condition 

d 3 A 3 l (x) = 0. (3.17) 

Since d 3 A' 3 = d' 3 A' 3 n this gauge implies A' 3 (x) = and hence leaves from A 3 (x) only the 
a^-independent part a 3 (x). 



By construction (see eqs. fl3.16Q and (|3.11|) ) the reduced Abelian fields a 3 (x) and e 3 



x 



are canonically conjugated variables. Note also, that the change of variables from A 3 

14 



to (A' 3 n ,a 3 ) (and correspondingly from E 3 to {E'^,ez)) does not yield any non-trivial 
Jacobian since A' 3 and a 3 are orthogonal coordinates in the sense that they belong to 
orthogonal subspaces of the Hilbert space of periodic eigenfunctions of id$. 



The gauge condition ( 3.17]) has also the advantage that it enormously simplifies the op- 
erator (|3.9|) 

Df = S ah d 3 + df (x) =: df (3.18) 

which enters the charged part of Gaufi law ( |3.4j ) 

r s = dfE\ + T\. (3.19) 

(The evaluation of the eigenvalues, and hence the inversion of d 3 becomes trivial since 
a 3 (x) is independent of x 3 , see below.) Let us also mention that equation ( |3.8|) and (|3.17|) 



define a variant of 't Hooft's maximal Abelian gauge fig] , which preserves invariance under 
^-independent Abelian gauge transformation. 

For the time being, let us assume that <i 3 6 has no zero eigenvalue in the charged subspace. 
(We will later see that the system dynamically avoids configurations a 3 (x) = giving rise 
to zero modes of d 3 .) The charged part of the Gaufi' law can now be used to eliminate 
the charged part of E 3 by using 

Hn = ^6 (e$ + (d^f rl) . (3.20) 

Later we will observe that the corresponding functional determinant det df' will be can- 
celled by the Faddeev- Popov determinant. 



We can use now the two constraints ( ft.15 ) and (|3.20| ) arising from the Gaufi' law to 



integrate out explicitly the electric field variables E' 3 C °,E^ leaving from E 3 only the x 3 - 
independent neutral part e 3 ( |3.11| ). Furthermore the two gauge constraints ( |3.8|) and 



( [3.17] ) eliminate the gauge variables Af 1 and A' 3 n respectively, leaving from the gauge 
potential A 3 only the neutral ^-independent part a 3 (x). Since the changes of variables 
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from E 3 ° to E 3 °, e 3 and analogously from A^ to A 3 °, 03 are trivial, i.e. the corresponding 
jacobians equal one. We then obtain from ( |2.30| ) 

Z\J\ = f V(A ± ,a 3 ,E ± ,e 3 )]j5(r(A))5( 1 l)]jDetM ab -[Det(d' 3 )Detdf}~ 1 

x x° 

exp |^ L J d 3 x (e 3 d a 3 - ^e 3 (x) e 3 (z)) + J d 4 xE a ± d A a ± 

-\jd 4 x \(dz\ch)Tf) 2 + (d'^Y'ff + E ± E ± + E 2 ] } . (3.21) 



Here 8 [f c °(A)j denotes the gauge condition necessary to fix the residual invariance under 
^-independent Abelian gauge transformations, which is left by the constraints ( |3.8| ) and 
( p.lTp . This residual gauge will be fixed in the following section when we resolve the 
residual GauB law j± — 0. 



IV. RESOLUTION OF THE RESIDUAL GAUfi' LAW 



The residual GauB law constraint 

can be used to remove the a^-independent part of E™ which is longitudinal in the 1-2-plane 
(i = 1,2) defined by 

e± := Vx^T-^ J dx 3 V ± El = IEI , (4.2) 

where A' ± is the two dimensional laplacian, Vj_Vj_, in the Hilbert space of periodic 
functions with the zero mode omitted. Its inverse is defined in the space of periodic 
functions by the Green's function 

(*'_[, *1) = (zll-^H*') = £ 1 m -K-<)¥ , n± = (ni , n2 ) (4.3) 

-1 \ ' n±^0 -1 
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which obviously satisfies periodic boundary conditions. Note that the longitudinal pro- 
jector I defined by eq. ( |4.2| ) is in fact an orthogonal projector, 1-1 = 1. This follows from 
the relation 

j d 2 xl(x\A ± \x")(x"\^-\x') = 6® ( x± - X 'J - J_ , (4.4) 

where 

k ± 

is the 2-dimensional periodic 5-function (x± + eiLi,y±) = 5^ (x±,y±)) and the last 
term in ( |4.4|) arises from the fact that in A' ± ( [4.3| ) the zero mode n\ = n 2 = is excluded. 
This term, however, does not contribute when I acts on vector fields Vi(x) periodic in x\ 
and x<i- In fact, from the definition of the longitudinal field e± (12) we find by using (|4.4|) 



V ± e ± = jj dx 3 V L E n ± - J d 3 xV ± E n ± , (4.6) 

where the last term vanishes for periodic electric fields, so that we obtain 

V ± e ± = 2j dx3 V ± E n ± . (4.7) 
The residual Gau8' law (|4.1|) then simplifies to 

1± = V ± e ± - p (2) = , (4.8) 

where 

-jJdx 3 ([A ± ,E ± ] n + 4) . (4.9) 



P (2) 



Since by definition, e± ( (4.2|) is a curl-free, 2-dim. vector field, Vj_ x e± = 0, it has the 
representation 

e± = -V ± <p{x ± ) , (4.10) 

where the scalar potential <£>(x) follows from the residual Gaufi' law 
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¥> (x±) = / d 2 y L G^ (x ± , y ± ) ( y± ) . (4.11) 



In fact, inserting ( |4.11 ) into (|4.10|) and taking the divergence we find with the help of 

(El 

V ± e ± = p^-p , p = -J- / rfW 2) = TT-rW ' ^ 412 ) 

where Q n is the total charge (in the Cartan subgroup), which according to ( |3.3|) has to 
vanish for periodic .Ef -fields, so that p = and (|4.10| ) solves, in fact, ( |4.8| ). 

For non-vanishing total charge Q n ^ Gau8' law requires one to abandon the periodic 
boundary condition to the electric fields and the second term in ( f4.6| ) no longer vanishes. 



Even in this case eq. ( |4.8| ) is still solved by eqs. Q4.10D , ( |4. 1 1| ) . Therefore the resolution 
of the residual part of GauB' law leads to the elimination of the longitudinal part e± of 
the neutral vector field E\ and we are left with the transversal part 

E' ± = Ej_-e± (4.13) 

as dynamical quantity. 



Note, that only the charged parts A c ^ and Ef 1 enter p^ ( |4.9[ ) and thus e±. Furthermore, 
e± and E'™ live in orthogonal subspaces of the Hilbert space of periodic functions in x 3 . 
Therefore the change of variables from {E 1 ^, E c ^j to (E 1 ^ = E'™ + e±, E c ^j does not give 
rise to any non-trivial Jacobian. Let us also emphasize that after resolution of GauB law 
( [4.8D e± is not an integration variable but a function of Ef 1 , A c ^ and independent of the 
remaining integration variables E'\ , A\ etc. 

We can exploit now the residual invariance under ^-independent Abelian gauge transfor- 
mations, left by the constraints ( |3.8|) and (13.171 ), to remove the field 



a ± (x) = (lA n ± ) (x) (4.14) 



canonically conjugated to e±. Since by definition of the longitudinal projector I ( [4.2|) this 
field is curl-free, x a± = and, for periodic fields A±(x), satisfies the relation (c.f. eq. 

(EZP) 
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V ± a ± (x) = ^Jdx 3 V ± A n ± (4.15) 
it suffices to require the gauge 

jJdx 3 V ± Al(x,x 3 ) = (4.16) 

to make a± vanishing 

a ± = . (4.17) 
In the following we will denote by A' ± the field satisfying the gauge condition (|4.16|) , i.e. 



A' ± = A±-a±. (4.18) 



Since e± and a± live in the Cartan subgroup we can trivially extend eqs. (|4.13|) and (|4.18|) 
to the charged field components, where they read 

E >ch = E ch ^ A'f = Af. (4.19) 
We can then express e± defined by ( |4.1U| ) and ( f4.11| ) as 



e± (*) = Vlj J d 3 yG^ (x ± ,y ± ) ([A' x , E' ± ) n + J ") (x°,y) . (4.20) 
(Note that only the charge fields Ef 1 , A c ^ enter the commutator). 

Similarly we can express the Gau8 law generator I\ (|3.4[) in terms of the new variables 



n = V± {E'l + e ± ) + [A' ± , E' ± } n + 4 , 

Tf = V±Ef + [A' ± , E' ± + ej + Jt • (4.21) 

Furthermore since e± (x) ( |4.20| ) is independent of x 3 , i.e. d 3 e± = d' 3 e 3 = 0, it drops out 
from 

fff'ri = ±%rt . (4.22) 

and the neutral part of the Gau8 law generator ( ^.21| ) can be replaced by 
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rj = V±E'l° + [A' ± , E' ± } co + J C ° 



(4.23) 



Recall, that the change of integration variables from E± — > {E'™, e±) yields a trivial 
Jacobian equal to one since E'™ and e± are orthogonal components of E\ in the Hilbert 
space of periodic functions in x 3 . The same is true for the change of variables from A\ to 
(A'™,a±). Therefore, after complete resolution of Gau8' law and implementation of the 
gauge fixing contraints, we are left with the following functional integral representation 
of Yang-Mills theory 

Z = I V (E' x , e 3 , A' x , o 3 ) JJ DetM ab Der l (d f , 



L J d 3 x (e 3 d a 3 - ie 3 e 3 ) + J d 4 xE' ± d Q A' ± + J 



a 3 J 3 + A\J' 



d' yA +e'e' 



(B(A')f 



(4.24) 



where the residual abelian gauge constraint (4.16) has been used to replace the perpen- 



dicular field A± by its two-dimensional-transversal part A' ± , see eq. ( 4.18 ). Furthermore, 
the magnetic field B(A') is defined in terms of the reduced field variables due to the 
implementation of Gaufi' law by 



,3 



LI 



D[,d 3 



d 3 = d 3 + a 3 , 4 = d' 3 + a 3 
D' = d t + A' 



(4.25) 



Let us also emphasize, that there are no cross terms between the reduced electric field 
E' ± and the static electric field e±. This is a consequence of / dx 3 V±E' ± = 0, which holds 
due to the periodicity of the fields. 



It remains to calculate the Faddeev-Popov determinant which is done in the next section. 
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V. EVALUATION OF THE FADDEEV-POPOV DETERMINANT 



For the above chosen gauge the Faddeev- Popov determinant is straight forwardly evalu- 
ated. The two abelian gauge fixing conditions ([3 .17 ) and (|4.16 ) are independent of each 
other, i.e. they belong to orthogonal subspaces of the Hilbert space of periodic functions 
in x 3 G [0, L\. Both conditions can therefore be absorbed into a single gauge constraint 
for the neutral component of the gauge field 

1 ) 

f a °(x) = 3 Af (x,x 3 ) + Vx-y dx 3 A a f (x,x 3 ) . (5.1) 



Furthermore the gauge fl3.8|) defines a color charged gauge functional 



f* = A*. 



(5.2) 



For the above gauge functionals (|5.1|) and ( |5.2| ) the Faddeev-Popov kernel A4 ab (x, y) 
becomes (x = yo) 



M ab "(x,y) = Df°(x)V y 3 5 {3) (x - y) + yDf°(x)V y ± 5 {2) (x ± - y ± ) 



L 



M a \x,y) = bf{x)5^{x-y) . 



(5.3) 



This expressions hold so far for arbitrary gauge field configurations. We need, however, 
these expressions only on the gauge manifold, i.e. for those field configurations which 



fulfill the above chosen gauge constraints. Using f ab ° c ° = 0, which implies 
Faddeev-Popov kernel reduces at the chosen gauge orbits to 



the 



M ab {x,y) 
( 
V 



V 



(5.4) 







M abo M ab 

_ 5 a b (yx V x 5 (3) (f _ + 1 V1V15 (2) [:r ± - u ± ; 

-±Af°{x)V x ± 5W ( x± -y ± ) df{x)5^{x-y) ) 

Since this matrix has triangle form, we find for the Faddeev-Popov determinant finally 
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DetM ab (x,y) = Det 



_ S a 0b0 ( V x V x j(3 ) + IviV15 (2) (x ± - y ± ] 



Det (dfS® (x-y)) . 



(5.5) 



It factorizes into contributions arising from the Cartan subgroup (first factor) and the 
coset space. The former one is an irrevelant constant and will be dropped in the following. 
The contribution from the coset space can be easily calculated since the eigenvalues of df 1 
are analytically known, see Appendix A. But for the moment we do not need the explicit 
form of Det d^. 

A glance at eq. (|5.5| ) shows that (the non-trivial part of) the Faddeev-Popov determinant 
cancels precisely the determinant (^Detd^j arising from the resolution of Gaufi' law. 
Consequently eq. (^4.24|) reduces to 



Z = j V {E' x , e 3 , A' ± , a 3 ) exp j-^ L J d 3 x (e 3 d a 3 - ^e 3 e 3 ^ 

+ J d 4 xE' ± d A' ± ~\jd 4 x [(4 -1 r±) 2 + E' ± E' ± + e 2 ± + (B'(A)) 2 ^ } . (5.6) 

This is the desired functional integral representation of Yang-Mills theory in uncon- 
strained, gauge-fixed variables, resulting from a complete resolution of GauB law. Note 
that in the unconstrained theory the functional integration over the canonical variables is 
performed with flat integration measure. (There is no preexponential factor, e.g. a func- 
tional determinant, which could be interpreted as non-trivial measure.) This is obviously 
a general feature of Yang-Mills theory in unconstrained variables (provided one chooses a 
gauge condition which is canonically conjugated to the GauB law constraint) and could, 
perhaps, have been anticipated in view of the fact that the Faddeev-Popov kernel is given 
by A4 ab (x,y) = [f a (x), r 6 (?/)}, where {, } denotes the Poisson bracket. 

The cancelation of the Faddeev-Popov determinant against the determinant arising from 



the resolution of GauB' law was also obtained in ref. 0], where the gauge was fixed by 
demanding that the antisymmetric part of the matrix E? vanishes. In that case GauB' law 
requires the vanishing of the antisymmetric part of A a { and one ends up with a functional 
integral over the symmetric parts of Af and Ef, where unfortunately the remaining electric 
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field variables cannot explicitly been integrated out. In this respect the present approach 



has the advantage over refs. [0,0 in that the remaining unconstrained electric field 
variables in ( |5.6| ) can still be integrated out in closed form. This will be done in section 
7. 

Before concluding this section let us notice that assuming a flat integration measure the 
functional integral representation (|5.6|) could have also been derived by starting from 



the Yang-Mills Hamilton operator in unconstrained variables obtained in ref. |30| in the 
canonical operator approach and following the standard procedure |37|]. In this sense 
the present functional integral derivation of the unconstrained Yang-Mills theory (|5.6|) is 



equivalent to the canonical operator approach of ref. ||30|| . We believe, however, that the 
functional integral representation derived in the present paper (|5.6|) is more flexible than 
the operator approach when it comes to an approximate solution of the theory. 

Finally a comment on the gauge-fixing is in order. We have fixed the gauge in such a way 
to remove the components of the gauge field Ai(x) which are canonically conjugate to 
those components of the electric field Ei(x) which are eliminated by Gaufi law. This has 
lead to the gauge conditions (|3~8| ), (|4.16|) and ( [3.17|) , which eliminate A^ h ,a± and make 



A3 = tig independent of X3. These gauge constraints do, however, not yet fix the gauge 
completely but leave a residual gauge invariance which consist of i) (global) permutations 
of the color indices of the fundamental representations, i.e. elements of the Weyl (sub-) 
group Sn of the gauge group SU(N), ii) global Abelian gauge transformations and iii) 
displacement transformations Q = e~ ax with a an arbitrary but fixed c-number 3-vector. 
These residual gauge symmetries were also found in ref. |JIJ. For completeness we work 
out the emergence of these residual gauge symmetries in the present functional integral 
approach in Appendix C. 
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VI. THE WILSON LOOP 



Below we evaluate the potential between two static color charges or, equivalently, the 
Wilson loop 

W(C) = (TrP exp (- j> dx^A^x)) ) (6.1) 



in the (Euclidean version of) the gauge- fixed theory defined by ( |5.6| ) . For simplicity 
we consider a planar rectangular Wilson loop C, which by Lorentz invariance can be 
placed into the — 3 plane. One should note here, however, that the present approach 
( |5.6|) has not been formulated in a manifestly Lorentz covariant way, although all Green 
functions (calculated in the full theory) will respect Lorentz invariance. As a consequence 
the quality of approximations will depend in general on the chosen Lorentz frame. 

The present approach obviously singles out the 0— and 3— axis. (It starts from the 
Weyl gauge and eliminates most of the degrees of freedom of A 3 and E 3 by gauge-fixing 
and resolution of GauB' law, respectively.) We therefore expect that the Wilson loop is 
most efficiently evaluated when placed in the — 3-plane. Then the A' ± field will not 
explicitly enter the Wilson loop. Therefore we will ignore it together with its conjugate 
variable E' ± since we anyhow expect the dominant infrared behaviour to be governed by 
the a 3 (x),e 3 (x) fields. The generating functional of axial-gauge fixed Yang- Mills theory 
Q5.6| ) reduces then to 

Z[J] = = J V (a 3 , e 3 ) exp j L J d 3 x (e 3 d a 3 - ^e 3 e 3 ^ + J a 3 

-yH j °^k j ° + w 2 }}}- (6 - 2) 

where 

ef = e ± \ E±=0 = -Vj± / d 3 yG^ (x ± ,y ± ) J%{y) . (6.3) 

The last two terms describe the interaction between static charges Jo- The last term can 
be cast into the form 
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=f) 2 = _ / dx d 2 x L d 2 y ± J^ (x , x ± ) G (2) (x ± , y ± ) J C ° (x , y ± ) 



(6.4) 



where a partial integration has been performed and 



J = — I dx 3 J (x) . 



(6.5) 



This quantity obviously vanishes in the thermodynamic (infinite volume) limit L — > oo 
for any localized charge distribution Jq{x). To illustrate the meaning of this term let us 
consider two opposite Abelian charges (q, —q) separated by a distance R. If we place 
these two charges on a line parallel to the 3-axis, e.g. 



(6.6) 



then obviously Jo (x) = and this term does not contribute. But it does contribute when 
we place the charges in the x — y plane, e.g. 



J^(x) = q co 6(x 3 )6(x 2 ) 



(6.7) 



If we take, for simplicity, the thermodynamic limit Li, L2 — > 00 of (x±, y±) ( |4.3| ) 



G {2) (x±,y±) = In \x± - y ± \ 



(6.8) 



we receive from eq. ( |6.4j ), besides an infinite constant, a logarithmically increasing poten- 
tial. However, we do not expect that eq. QS.2p, which discards all perpendicular degrees 
of freedom A' ± , E' ± , can give a realistic description of the interaction between two charges 
sitting in the 1-2-plane. As discussed before the present approach singles out the 3-axis 
and in fact preserves the rotational symmetry around the 3-axis. Let us therefore consider 
the axial symmetric charge distribution (6.6). In this case = 0, and from the second 
to last term in ( |6.2|) we obtain the static interaction potential 



(6.9) 



where G^ 1 ' is the Green's function of — d' 2 . If we again take the thermodynamic limit 
L — > 00 
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(x 3 ,y 3 ) = (x 3 \-^\y 3 ) = l -\x 3 - y 3 \ (6.10) 
we obtain a linearly raising potential 

V = aR (6.11) 

with a string tension 

° = j 2 m)? ■ (6.12) 

Here it is understood that 5(0) is regularized in an appropriate way. The above obtained 
interaction potential is in agreement with the findings of the canonical quantization ap- 



proach pOfl , see also ref. pi ] 



The emergence of the linear confinement potential in the 3-direction should come as no 
surprise since, except for the dummy x±, ^-dependence, eq. ( |6.2| ) represents the gener- 
ating functional for 1 + 1 dimensional Yang-Mills theory, which is known to confine. In 
fact, if we ignore the term (which, as seen above, vanishes in the thermodynamic limit 
L 3 — > oo for any localized charge distribution Jq(x) and furthermore depends only on the 
"dummy" coordinates X\,X2) and linearize the term quadratic in Jq by means of a field 
ao(x), and furthermore perform the integration over e 3 we obtain 

Z[J] = J V (oq, a 3 ) exp |^ J d 2 x± J dx dx 3 Q [d a 3 ) 2 

+ ^(d 3 a o y + a J + a 3 J 3 ^J |. (6.13) 

Here the first two terms in the bracket combine to f? p , where 

fp,p = djxdo - dpd/i , p, v = 0, 3 (6.14) 

and additionally a 3 satisfies, by construction (see. eqs. (|3.8|) , ( |3.16|) ), the gauge 

af = , d 3 a^ = 0. (6.15) 

In D — 1 + 1 the corresponding Faddeev-Popov determinant is an irrelevant constant. 



Thus eq. ( 6.13 ) represents in fact the properly gauge- fixed generating functional of 2- 



dimensional Yang- Mills theory, except for the parametric Xi, x 2 - dependence of the fields. 
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Eq. ( |6.13| ) can be regarded as the strong coupling limit of Yang-Mills theory. This 



interpretation is consistent with the result of ref. where a strong coupling expansion 
of Yang-Mills theory was performed and the leading order was found to be given by D = 2 
Yang-Mills theory. This result is also confirmed in the field strength approach |[38||. 



It is now straightforward to evaluate in the reduced, 2-dimensional Yang-Mills theory 
( |6.13| ) a Wilson loop in the 0-3 plane. One finds the area law in agreement with the linear 
rising potential between static charges as found above. 



VII. ELIMINATION OF THE ELECTRIC FIELDS 



In the gauged fixed Yang-Mills theory, where the Gaufi' law constraint has been fully 
resolved, the electric field variables occur still only quadratically in the exponent, so that 
these variables can be integrated out. The integral over e% is trivial. To perform the 
integral over E'j_ it is convenient to introduce a more compact notation. We define the 



kernel^ 



K ab (x,y) 



K~ al 
Furthermore, we define 



V 



_d>2fia b 





(44) 



\ 



ab 



6 {4 \x-y) 



J 



(7.1) 



(E' ± (x) + e ± (x)) c = I d i yP cc '(x,y)Ei(y) + e 



,(0)c 



where is defined by eq. (|6.3|). and 

if' = 6 cc '6 t j6^\x, y) - <f c » (v~) \ J dx 3 Af(x)6 bc '5^(x - y) 



(7.2) 



(7.3) 



This quantity fulfills the relations 



5 Note that, up to an irrelevant constant, Det K gives the square of the Faddeev-Popov deter- 
minant (|T^) . 
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V^(x,y) = V J f b 6^(x-y) 
Vj^'foy) = jf dx 3 D'l° b °(x)6W(x,y) . 



In this notation we have 
d 4 x 



d x 



r i K^Y 



J d A x [(D' ± (E' ± + e±) + p) K~ l (p' ± (E' ± + e ± ) + 
J d 4 x (D' ± (PE' ± + ef) + J ) K- 1 (D' ± (PE' L + ef ) + J c 
/ d'x ((E' ± P T + ef) (-D' ± ) + J ) (D> ± (PE' ± + ef ) + J c 



Here, we have also introduced the transposed projector P T by 



(PE' ± nx) = (E f ± P 1 ) (x). 



In this notation we can also write 



/ (E' ± E> ± + e\) = J (E' ± + e ± ) 2 = J (PE' ± + 



i 0) )' 



J (e' ± P t PE> ± + 2ef ) PE' L + (ef ) : 



(7.4) 



(7.5) 



(7.6) 



(7.7) 



In eq. (|7.5|) and below it is understood, that in D' ± = Vj_ + A' ± the 2-dimensional gradient 
operator Vj_ is replaced by the corresponding operator Vj_ with the zero mode excluded 



(x ± \V' ± \x') = i J2 n^M*-*') 



n± = (ni,n 2 . 



{7.1 



This is admissible since V±E± = V' ± E±. The integral over E' ± can then easily be carried 
out yielding for the transition amplitude^] 



6 Note in eq. ( |7.9| ) from D'j_e^ = V±e ® + A' ± e ® the term Vj_e_^ can be dropped since this 
quantity does not depend on xj, and hence vanishes when acted on with K^ 1 ~ (c^ 2 ) 1 
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L^J (f'.r {<) U (I {.r)f I. I d 3 xa 4 (.r)j: l (.r) 



Z[J] = [ V(A' ± ,a 3 )~[[Det- 1/ mexp\- 
+~ J d 4 x (ef - d A' ± + (D'J? + J ) K~ l D ± ) PR-'P T (ef - d Q A' A 



-D'K' 1 (D^ + Jo 



+ / d*x (efd A' ± - \ (ef ) 2 ) - \ J d*x (LrJ? + J ) i (i^ef + J 
-i/rf 4 x(5(A')) 2 



(7.9) 



Here, we have introduced the kernel 



H = P T P + P T f —D\ 4tD\ ) P 



K 



(7.10) 



which is defined in the space of the spatially periodic 2- dimensional vector functions. 



It is instructive to consider the limiting case A' ± = 0, where the kernel ( |7.10| ) reduces to 



H = 1 - V ± K~ 1 V ± . 

In this case the terms of the action ( |7.9| ) containing the static source Jo reduce to 

-1 



(7.11) 



Jo^D' ± PH-"P T (-D' ± ) 1 Jo - J ^ Jo 



J 



0" 



A' ± =o 'K-A' A 
and in the limit 03 = we recover the familiar Coulomb law 



Jn 



(7.12) 



J 



K — A' 



-J, 



Jo^Jo , 



A' = Al + d 2 3 . 



(7.13) 

^-L J a 3 =0 

Note that in the continuum limit the zero eigenvalue of A' disappears and A' — > A. In 
the Yang-Mills theory the infrared singular behaviour of the static Coulomb law is avoided 
by the presence of the abelian field 0,3(0;). The infrared behaviour of Yang-Mills theory, 
and in particular the confinement mechanism, should therefore be essentially determined 
by this fluctuating field. This is in agreement with the findings of the previous section. 

In Appendix B it is shown that for A' ± = 



Det' 1/2 H = Det 1/2 K ■ Der lh (K - A' ± 



(7.14) 



Here 
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Det 1/2 K = const. Det 1/2 (-44) = const • J (7.15) 

is precisely the Faddeev-Popov determinant ( |5.5|) which, as shown in Appendix A, up to 
an irrelevant constant coincides the Haar measure of SU(N) 
j rr 1 2 a k (x)-ai(x) 

J= [[j^ SmL o ,2^OC k {x)=0, ( 7 - 16 ) 

k>l L Z k=l 

where afc(x) are the diagonal elements of i a 3 (x) (c.f. also ref. |3"5|]). 



Note, that after imposing GauB' law we have obtained here the Haar measure of the gauge 
group for the functional integral over the gauge field as(x) although we started from the 
usual functional integral representation with flat integration measure but gauge fixed by 
the Faddeev-Popov method. This result is a manifestation of the observation [33] that 



the standard functional integral repesentation of Yang-Mills theory with gauge fixed by 
the Faddeev-Popov method fully respects gauge invariance even in the non-perturbative 
regime and that, in particular, the Haar measure of the gauge group (necessary for a pro- 
jection onto gauge invariant states) naturally arises from the Faddeev-Popov determinant. 



The Haar measure and hence Det 1/2 H vanish for degenerate field configurations 03 (x), 
for which two diagonal elements coincide, i.e. (03 {x)) kk = (a 3 (x)) u for k ^ I. 



Det' 1/2 H = exp ( — Tr log H) (7.17) 



Since 



2 

we receive an additional contribution to the effective action of the remaining (physical) 
degrees of freedom, which represents an action barrier to keep the system out of the 
singular field configurations. Such (energy) barriers have been also found in alternative 



formulations of gauge theory in terms of gauge invariant variables P|,p5[ • 



Finally, let us make a few comments concerning the relation of the present approach with 
those interpreting confinement as a dual Meissner effect arising from monopole conden- 
sation. In fact, the gauge defined by eq. ( |3.8|) and ( 3.17] ) is a variant of maximal Abelian 
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gauge [T^JT^]. In these gauges monopoles arises at those singular points in configuration 



space where the gauge fixing is not unique. The gauge ( |3.8| ) is not unique at those sin- 
gular points x = xs, where the field as(x) is degenerate, i.e. two eigenvalues of a^x) 
coincide. It is straightforward to show |T5| , |I6|| that near these singular points the gauge 
transformation Q (x) necessary to fulfill the gauge ( |3.8| ) is such that the Abelian part of 
Qd^Q^x) developes a "magnetic" monopole in fi = 0, 1,2 space. Field configurations for 
which the gauge fixing is not unique give rise to zeros of the Faddeev- Popov determinant. 
In fact in the present case the Faddeev-Popov determinant 

DetiD 3 ~ J ~ VDetK (7.18) 

vanishes at the singular points of degenerate 03 (x) field configurations. The field configu- 
rations of vanishing Faddeev-Popov determinant define the Gribov Horizon, which in the 
present context is therefore built up from monopoles. 



VIII. CONCLUDING REMARKS 



In this paper we have considered D = 3 + 1 dimensional Yang-Mills theory defined on 
a spatial torus. Using a variant of 't Hooft's maximum Abelian gauge (see eqs. ( |3.8| ), 
( p.!7| ), ( J4.16 )) we have performed a complete resolution of Gaufi' law. This has resulted in 



a functional integral representation of Yang- Mills theory which is entirely defined in terms 
of unconstrained, gauge fixed variables. These are the spatial gauge fields i — 1,2 
defined by eqs. ( (4.2|) , (^4.14|) , (|4.18|) , the neutral x 3 -independent part of which is transverse 
in 1-2-plane (see eq. (^4.16|) ). In addition an Abelian ^-independent field a 3 (x), defined 
by Q3.16Q arises, so that the total number of degrees of freedom is that of 2 (A 2 — 1) fields. 



This is the correct number of unconstrained degrees of freedom of massless Yang-Mills 
theory with gauge group SU(N). In this respect the present approach is more efficient 
that the pioneering works of refs. |19,20| where twice as much integration variables remain, 



since in that case the electric field variables cannot be integrated out in closed form. 
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We have worked here in the Hamilton functional integral formulation which obviously 
violates Lorentz covariance. Furthermore, the adopted gauge also violates spatial 50(3) 
invariance but preserves axial symmetry. In this respect the present gauge is advantageous 
over the gauge used in ref. jllj which violates also axial symmetry. Of course a complete 



elimination of all gauge degrees of freedom without violating any space-time symmetry 
would be preferable. This has been partly achieved in refs. PHEffil for 577(2) in D = 3, 4. 
This appraoch works in the canonical Hamilton (operator) approach, which obviously 
violates Lorentz covariance but preserves all spatial symmetries. Unfortunately there is 
no direct way to extend this approach to higher gauge groups SU(N > 2), although some 
attempts have been undertaken for SU(3). For SU(2) in D = 3 a complete covariant, 
gauge invariant description has been achieved in the so-called field strength approach [[E 



at the expense of a doubling of the degrees of freedom p8| . We consider, however, the 
violation of a global symmetry, which can easily be restored, a minor problem. Of course, 
the exact Green functions preserve all space-time symmetries even in gauges which violate 
these symmetries. 

As was illustrated in sect. 6 the field a^{x) represents the dominant infrared degrees of 
freedom, which in particular are responsible for the emergence of the area law. As a first 
step one might include only this Abelian field for a study of the infrared sector of QCD. 
In fact, due to our adopted gauge this is in the spirit of the Abelian dominance observed 
in lattice calculations performed in maximum Abelian type of gauge [ED . 
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APPENDIX A: EVALUATION OF THE FADDEEV-POPOV DETERMINANT 

Below, we evaluate the functional determinant of id 3 . Consider the eigenvalue equation 

iD'f<p b v = , (Al) 



where (p v has to satisfy periodic boundary conditions in x 3 . Multiplying this equation 
with the generators in the fundamental representation T a and using 

T a Df V b =[D,^} , ( p = (p T a (A2) 

the eigenvalue equation becomes 

% M + [»3^ (l/) ] =/W H - (A3) 
This equation is easily solved, since ia 3 = ia^ T C0 is a (traceless) hermitian diagonal 

N 

matrix, with real elements a k = (a 3 ) kk satisfying J2 a k = 0. Hence, the eigenvalue 

k=i 

equation reads explicitly 

id^ti + («* - a i) = »uVu • (A4) 
The periodic eigenfunctions are given by (v = (c,n)) 

= rfc***™ , U - = ^T ' ^ = °' ±1 > ±2 --- > (A5) 

where the r) kl denote the vectors of the Weyl basis of SU (N) , which with c—(r,s), r,s = 
1, 2, . . . , N is defined by 

»7« = ^<W* S • (A6) 
The corresponding eigenvalues read 

fJ> v = Hc,n = u; n + a r - a s c= (r, s) . (A7) 
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For fixed n, there are 2 eigenvectors <^ = ( r > s ) corresponding to the off-diagonal elements 
r ^ s. The corresponding eigenvalues come in pairs iv n ±\a r — a s \. Since trip^ = there 
are only N — 1 independent eigenvectors ip u=( - r ' s ^ with r = s, being degenerate with 
eigenvalue cu n . The total number of independent eigenvalues (for fixed N) is of course 
N 2 -l. 



We therefore obtain for determinant under consideration 

oo 

Det (iD'^j = Yl fi u = Yl Y[ /x c o in 



c n=— oo 



n n (^n + «r - « S ) 



= Det (id^) ■ Det (d 3 ) . (A8) 

Note that by definition of Z)' 3 the mode n = has to be excluded for r = s, which is 
indicated by the prime. The expression in the first bracket yields an irrelevant (diverging) 
constant, which can be absorbed into the renormalization of the functional integral. 



Using 



oo / 

sin a; = x Y[ I 1 

n=l V 



X 



(A9) 



the expression in the second bracket can be transformed to 

oo 

Detidz = Y[ ( a r — oi 3 ) J| [(a r — a s ) 2 — to. 

r^s n=l 



const' Y\ («r — oi s ) II 1 — ( 

r^s n=l \ 



a r — a s 



= const 



const Y\ 



n2 . cx r ot s 



2 . a r - a s 



r>s 



(A10) 



Thus up to an irrelevant constant this determinant agrees for Y.r=i a r — with the Haar 
measure of the group SU(N) 

(a k - ai) 



J {La) = Yl sin 2 L- 



(All) 



k>l 
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APPENDIX B: EVALUATION OF DET H 



In what follows, we work out the functional determinant Det H for A' ± = 0. In this limit 
H is block diagonal in color space, i.e. it has no matrix elements between the neutral and 
charged color space 

! fjaok) Q \ 



H 



ab 



Bl 



i H at j 



This is because the same is true for the matrix K ( |7.1|) . Therefore, the functional deter- 
minant Det H factorizes as 

Det H = Det {n) HDet {ch) H . (B2) 
Since, the color neutral part of H is given by 



H?f° = 5 aobo ( % + Vj-L V; ) . (B3) 



Det^H is an irrelevant constant, which will be ignored in the following. 

For the evaluation of Det^ ch 'H we consider the corresponding eigenvalue equation 

Hp^x) = im%{x) , (B4) 
which with the explicit form of H ab reads 

- VJ (R- 1 )"" V l 4>%x) = 0* - !)#(*) . (B5) 



Here, the eigenf unctions have to satisfy periodic boundary conditions and we have intro- 
duced the abbreviation K = —^3^3, which is the matrix K ( |7. 1| ) in the charged subspace. 



The eigenfunctions <j)\{x) represent 2-dim. spatial vectors, which we can split in longitu- 
dinal and transverse parts 
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satisfying 

d l rt(xf = , d^xf = d i( f>l(x) . (B7) 

Obviously any (spatially) transverse vector function (t)\{x) T gives rise to a zero eigenvalue 
of Vjif _1 Vj and hence to an eigenvalue \i = 1 of H. 

The longitudinal part (p^(x) L gives rise to a non-trivial eigenvalue |i / 1. For these 
eigenvalues eq. (|B5|) can be simplified. 

Operating on equation (|B5|) from the left with V£ and defining 

V \x) = V'M{x) , (B8) 
the eigenvalue equation becomes 

-A' ± (K- i y b <f?=fr-l) V P. (B9) 
Therefore the non-trivial part of the determinant of H is given by 

DetH = Det (l - A'^fsT 1 ) . (BIO) 

Note that the kernel of the r.h.s is a matrix in color and functional space but a scalar 
in ordinary space, contrary to H which is also a matrix in the 2-dim. Euclidean space 
spanned by the £j =1 2 -axis. The missing dimension on the r.h.s. is due to eigenvalues 
fi = 1 of H . For later use it will be convenient to separate off K~ x yielding 

Det (R - A', ) 

DetH = DetK ■ (BU) 
Since id^ is the hermitian operator K = —d\ is positive semi-definite, while (— Aj_) is 
strictly positive definite. Therefore, (^K — A' ± ) is a positive definite operator and con- 
sequently Det~^ — A'j^j is non-singular, even for field configurations a 3 (x) for which 
K has zero eigenvalues. These field configurations do not contribute to the transition 
amplitude due to the presence of Det^ K. Note, that this determinant agrees with the 
Haar measure of SU(N) given in Appendix A. 
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APPENDIX C: RESIDUAL GAUGE INVARIANCE 



The spatially periodic boundary conditions to the gauge fields Ai(x) restrict the gauge 
transformations Ai — > Af to those gauge functions Q(x), which satisfy the equation 



A, ^k(x 



0, (CI) 



where 

Cl k (x) = fi f (x + Le k )tt(x) . (C2) 

For simplicity we have set here L\ = L2 = A = L. This equation, which has to be 
fulfilled for all Ai(x) and all k = 1, 2, 3, is solved for 

n k ( x ) = z nk , (C3) 

with 

Z n = e'^ = e i2 7r ; n = 0, 1, 2, . . . , JV - 1 (C4) 



being an element of the center of the gauge group. The center of the group is defined by 
the set of elements commuting with all elements of the group. Here n is an element of 
the Cartan algebra TC. For SU(2) the center of the group is given by 

0„=o = 0,Z = 1 
<j) n =i = -inr 3 , Z x = -1 . (C5) 



Eqs. (£D, (|CJ imply 



Q (x + e k L) = Z*n(x) . (C6) 



The gauge conditions chosen above (see eqs. ( |3.8| ), (|3.17|) , ( |4.16| ) ) do not yet fix the gauge 
completely. There is a residual gauge symmetry left, which will be exhibited below. 
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First we note, that all three gauge fixing conditions are left unchanged under permutations 
of the color indices k, I of the fundamental representation for the gauge group SU(N). 
Such permutations of the basis color vectors are generated by those global gauge trans- 
formations S G SU(N), which are N-dimensional matrix representations of the symmetry 
group Sjy (group of permutations of N elements). 

These transformations form the Weyl group. For SU(2) the Weyl group consists of two 
elements Sq = 1 and Si = —it 1 = 2T 1 , which correspond, respectively, to the trivial 
permutation and to an exchange of the two color indices. The non-trivial permutation in 
fact represents a rotation in color space around the 1-axis through an angle ir 

S l = -it 1 = e~ 1 ^ = e wTl . (G7) 

Note, that Si is not an element of the Cartan algebra. These considerations can obviously 
be extended to larger gauge groups SU(N > 2). 



Besides the above discussed global residual gauge symmetry the gauge conditions (p.8j) , 
( |3.i71 ), ( |4.16|) are also invariant under the so-called displacement symmetry |K],[ftJ. Let 



us show how this comes about. The gauge condition Af 1 = leaves residual gauge 
transformations of the form 

Q = U(x)S , U{x) = e~^ x) , u{x) G H , (C8) 

where S G and U{x) is a gauge transformation in the maximal abelian subgroup 
(invariant torus). 

The second gauge condition d^A^ = is left invariant by gauge transformations of the 
form ( |C8j ) provided that 



d^u{x) = . (C9) 

This equation is satisfied if u{x) is of the form 

u(x) =cu (0) (x) +uj^(x)x 3 ; J°\x),u^\x)en. (CIO) 
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The quasiperiodic boundary condition to the gauge function (|C6|) requires 

This equation has to be fulfilled for all x and is solved for 

co W (x) = i(2vr^3-0„ 3 ) -«3 (C12) 

with k 3 = dia f/4 , fc^ , • • • ^3 ) being a N-dimensional traceless diagonal matrix, 

Y, kjp = 0, with integer entries k®. 

i=i 

Finally, the third gauge constraint restricts the residual gauge transformations to such 
function satisfying 

V±V±(u,n + ±Lu,l 1 \z))=Q, (C13) 



which in view of equation QU12| ) reduces to 

V±V ± u/°)(x) = 0, (C14) 

which implies 

cj( )(x) =/?(x°) + , f3,a ± eH. (C15) 

The quasiperiodic boundary condition (|C6| ) requires the Q!j(x ) to fulfill the relation 

e - a± L = ^ = j (C16) 

which is satisfied for 

a ± = j(2^k ± i-<p n± ) (C17) 

with = dia (k±\ k±\ . . . , fc^ ) being again a traceless diagonal matrix with integer 
entries , satisfying J2 = 0. 

The residual gauge symmetry left after the three gauge constraints have been implemented 
is therefore given by 
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Q( x ) = Se-W+s® . (C18) 

Besides the discret symmetry S generating the permutation of color indices the residual 
gauge symmetry consists of a global abelian gauge transformation e _/3 and the so-called 
displacement symmetry e~ ax . 
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